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Hydrodynamic surface waves propagating on a moving background flow 



•^ experience an eff^ective curved space-time. We discuss experiments with 

I 

^ gravity waves and capillary-gravity waves in w^hich ^ve study hydrody- 

• namic black/white-hole horizons and the possibility of penetrating across 

O 

•^ them. Such possibility of penetration is due to the interaction with an 

r-| additional "blue" horizon, which results from the inclusion of surface 



tension in the low^- frequency gravity-w^ave theory. This interaction leads 



^ to a dispersive cusp beyond which both horizons completely disappear. 

^^ We speculate the appearance of high-frequency "superluminal" correc- 

^»o tions to be a universal characteristic of analogue gravity systems, and 

^ discuss their relevance for the trans- Planckian problem. We also discuss 

I the role of Airy interference in hybridising the incoming w^aves with the 

<^ flowing background (the eff^ective spacetime) and blurring the position of 

V^ the black/white-hole horizon. 



I. INTRODUCTION 

Several physical systems reproduce certain properties of astrophysical objects like 
black holes: they exhibit an effective curved space-time when a wave propagates in a 
"moving" medium [I1I2]. Examples can be found in acoustics, dielectrics, optical fibres, 
micro- wave guides, Bose-Einstein condensates, superfluids, ion traps... Starting with 
the seminal works of White |3], Anderson & Spiegel |1], Moncrief [5] and Unruh [6j, 
there has been a growing interest for such analogue models of gravity in order to 
simulate and understand the physics of wave propagation on a curved space-time. 
The case of interface and surface waves was initiated by Schiitzhold & Unruh, who 
derived the equation of propagation of long gravity waves moving on a background 
flow in terms of a general relativistic metric [7j. A black hole, or its time-inverse: a 
white hole, can indeed be mimicked by the interaction between interfacial waves and a 
liquid current [7HIl|. Just like the expected behaviour of light near the event horizon 
of a black hole, gravity waves cannot escape a hydrodynamic black hole featuring a 
trapping line caused by the velocity gradient of a sufficiently strong background flow. 
In hydrodynamic experiments, it is usually more convenient to simulate a white hole. 
Then, gravity waves cannot enter a hydrodynamic white hole featuring a blocking line. 

In a previous series of experiments, such an artificial white hole was created and 
observed in a laboratory, using water waves in the presence of a counter-fiow inside 
a wave-tank with varying bottom profile [8]. Here, we provide a status report on a 
new series of ongoing experiments, in which we study the behaviour of surface waves 
of various frequencies in the gravity and gravity-capillary regime in the presence of 
such white-hole horizons. We encounter a good qualitative agreement with the theory 
developed in |9], |10] , and interpret the interesting quantitative differences. We are in 
particular interested in the following key points of the theory of gravity-capillary surface 
waves. When including surface tension, a second and third horizon appear on top of the 
white-hole horizon for gravity waves: a "blue" horizon associated to mode-converted 
blue-shifted waves, and a "negative" horizon. This negative horizon is associated to 
the appearance of negative-energy waves [8] (recently observed in a similar setup [TT]). 



an essential feature of Hawking radiation (the quantum glow of super-massive objects). 
The blue horizon can interact and even merge with the white horizon providing two 
scenarios in order to escape an artificial black hole. One scenario, consisting of a dou- 
ble bounce with mode conversion, was experimentally discovered nearly three decades 
ago by Badulin et al. [12] , and interpreted in terms of the black/ white- hole analogy in 
|10j . The other scenario is novel and consists of a direct dispersive penetration across 
the white-hole horizon. We have observed strong indications to validate this second 
scenario, and hope to confirm these quantitatively in the near future. This would 
establish that one can enter into the (normally forbidden) white-hole region by send- 
ing high-frequency capillary waves. These are not blocked at the primary white-hole 
horizon, contrarily to the low-frequency gravity waves, nor even at higher counter-flow 
velocities]^ Indeed, the horizon completely disappears above a certain frequency due 
to the merging and consequent disappearance of the blue and white horizons. 

These scenarios are not limited to hydrodynamic models, but are generic conse- 
quences of the dispersive properties beyond the relativistic regime. A double-bouncing 
scenario is possible in any system where "subluminal" dispersion (group velocity Cg 
decreasing with the wave number k) at intermediate wavenumbers k gives place to 
"superluminal" dispersion {cg increasing with k) at higher k. This latter condition 
(superluminal high-fc dispersion, irrespective of the behaviour at intermediate k) alone 
is sufficient for direct dispersive penetration. 

II. PRELIMINARIES 

Water waves in the presence of a uniform current are described by the dispersion 
relation [T2HT5] 

{u - U.k)2 = (gk + -kA tan\i{kh) (1) 

where u/27t is the frequency of the wave in the rest frame and k the wavenumber; g 

denotes the gravitational acceleration at the water surface, p the fluid density, 7 the 

^ In the presence of an ever-increasing counter-flow velocity, the capillary waves will continuously 
blueshift and ultimately vanish through viscous damping. We will come back to this point in 
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surface tension, U < the constant velocity of the background flow and h the water 
depth. For water, p = 1000kg. m"^ and 7 = 0.073N.m~^. The flow induces a Doppler 
shift of the pulsation u. 

The relativistic regime with its Schwarzschild-like metric corresponds to the shallow- 
water limit [kh ^ 1) of gravity- wave propagation [7]: [u — U.k)^ ^ ghk"^, with the 
relativistic "invariant" speed c = \^gh. An analogue black/white-hole horizon will then 
appear when |t/| = y/gh. But the concept of horizon can easily be generalized to any 
location where the group velocity Cg = duj/dk vanishes. The shallow-water limit kh <^1 
can be realized for example in the circular hydraulic jump, which spontaneously forms 
a hydrodynamic white hole [TB] . In wave-channel experiments, both the shallow- water 
and the deep-water limits can be probed, depending on the setup. Interaction with 
a counter-current will lead to a blueshifting of the incident waves: k increases. But 
variations of the counter-current flow rate are typically achieved by a variation of h 
(e.g. through the immersion of a bump). These two effects compete and various kh 
regimes can in principle be realised. At sufficiently high wavenumbers k, small-scale 
dispersive corrections appear due to capillarity, which will turn out to be crucial to 
cross a horizon. Since this is our main object of study, we will from now on focus on 
the deep-water limit {kh ^ 1), and start again from long- wavelength gravity waves, 
proceeding step by step towards higher k. 

When a gravity wave meets a counter-current, the incident wavelength diminishes 
and the wave height increases [T2HI11 [T7H^ . According to the ray theory, the wave 
amplitude would diverge when blocking occurs. However, such caustic for the energy 
is avoided by a regularization process. Due to the velocity gradient, the incoming 
waves are somewhat diffracted before being stopped at the blocking point, where Cg 
changes sign and the waves are reflected. Blue-shifted modes are therefore created 
through a process of mode conversion at the blocking point. Since the incident and 
blue-shifted waves have the same wavenumber at the blocking point, they interfere and 
a spatial resonance appears. The diffraction implies that the flgure of interference is 
not a simple standing wave. In the deep-water limit valid for our experiments, an Airy 
interferences pattern appears |9], fT9ti2T] . The blocking point itself is a saddle- node or 



tangent bifurcation [9J: it marks the point where the two real solutions disappear. 

When including surface tension [10], the deep-water {kh ^ 1) dispersion relation 
becomes 



{uj - Vky ^gk+ -^k-^, (2) 

or (tu — Uk)^ ~ (yf/c (1 + I'^k'^), where Ic = ^ /— is the capillary length (/^ — 2.7mm for 
water). The gravity waves are still blocked, with the blocking velocity for pure gravity 
waves (7 = 0) given by \Ug\ = gT /87c, with T = 27c/u) the period. Moreover, the 
blue-shifted waves are also stopped at a new blocking point on their backward drift 
(group velocity Cg < 0, along with the background flow U, although the phase velocity 
c^ = uj/k > 0). The capillary asymptotic limit for the blocking velocity of these blue- 
shifted waves is U^^^o = Uy = — v^ ( — j [lOJ. At this second blocking point, the 
blue-shifted wave merges with a new capillary solution, which appears (again through 
mode conversion, see Fig. IT]) at this secondary turning point. The capillary waves 
propagate in the same direction as the original incident gravity waves. These newly 
created capillary waves are not blocked by the primary saddle-node line but go through 
the gravity horizon [TO]. Badulin et al. observed experimentally that gravity waves can 
undergo such double bouncing behaviour followed by conversion to capillary waves 
which propagate into the forbidden region, and finally vanish by viscous damping [T^ . 
Fig. [1] illustrates the process. In the context of the black hole analogy, when time- 
reversing these observations, one concludes that incident long-wavelength gravity waves 
cannot escape from a trapping region (black hole) unless they are converted into the 
capillary range. 

The second escape route consists of creating incident waves from the start in the 
capillary range. The horizon completely disappears below some critical period T^, 
determined by the cusp formed through the interaction of the white and blue horizons. 
Waves with T < Tc can then propagate straight ahead, avoiding any horizons, and enter 
the "forbidden" region (or escape from the trapping region). The cusp can clearly be 
identified graphically from the {U* vs T) phase diagram for deep-water waves {kh ^ 1), 
see Fig. [2| where U* represents any critical or blocking velocity. 
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FIG. 1. Dispersion relation (pi) for varying counter-flow U (the curves rotate clockwise with 



increasing \U\), plotted in the form uj = Uk it \ igk + -k^j tanh{kh). The group veloc- 
ity Cg = -^ corresponds to the slope of the green/blue curves (green: positive co-moving 
frequency to' = u — U.k, blue: negative u'), and horizons are characterized by local min- 
ima/maxima. The double-bouncing observed by Badulin [12j corresponds to the dashed red 
line (see also [10]). 



The theoretical phase diagram Fig. pi was derived in [TO]. We briefly recall the main 
steps in its derivation. We start from the cubic dispersion relation (pi) for water waves: 
{u — UkY ^ gk + -k^, where we have taken k > 0. A double root k2 of this cubic 
equation, characteristic of a saddle-node point (and hence a horizon or turning point), 
is such that {k — ki){k — A;2)^ = 0, where ki is the remaining simple root. 

After some straightforward but tedious algebra, this constraint leads to a quintic 
equation for the critical velocity U*, corresponding to all possible resonances (for k > 
0): 



4pV[f/^ + ]^U^ + ^—U' - ^^f/^ - 6^t/ - 



4w 



P9 
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7^^ 277V 
pu 4 p'^g 



)] = (3) 



A dual quintic is obtained for /c < by reversing the velocity U -^ —U. 

Using this quintic, one can numerically compute and plot the velocity for saddle- 
node bifurcations or blocking points as a function of the period of the incident waves, 
as in Fig. pj Here we are mainly concerned with the gravity-wave blocking speed Ug 



and the blocking speed for blue-shifted waves (of which U^ is the asymptotic limit as 
T —7- oo), but we note that there also exists a blocking line for the negative-frequency 
waves [TO] . 

The well-known result Ug = —g/4uj in the gravity-wave limit immediately follows 
from eq. ^ when setting 7 = (or |f/| — t- 00). The asymptotic capillary limit 
U-y = —\/2i — ] is likewise obtained for T = ^ — )■ 00. From eq. ([3]), several 
analytic approximations can also be obtained. For example, to find an approximate 
expression for the blue horizon (the blocking of blue-shifted waves), we keep the terms 
up to first order in u and introduce the perturbative development 

[/-[/, (1 + ^) (4) 

in order to write: 



Solving for e, using a; — )■ 0, we obtain 



M 1 



- ^-^ (6) 



or 

f/* ^ f/, + ;^, /^ = -v^ (^-1) ^'^ + 27r|. (7) 



.^^—J— = -V2[—\ +271^. 
' TVgp \pj T 



Note that an identical computation for the negative quintic leads to a simple change 
of sign in the second term: to first order, the negative-frequency blocking line has the 
same departure (but in opposite direction) from the capillary asymptotic limit as the 
blue horizon (for T > T^). 

The cusp {Tc,Uc) = (0.425s, —0.178m.s^^) corresponds to the intersection of two 
saddle- node lines (the white and blue horizons), and is associated to an inflection point 
of the dispersion relation ([2]). In dynamical-systems theory, we anticipate a so-called 
pitchfork bifurcation [27]: the merging of two saddle- node bifurcations is equivalent 
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Ugl=gT/8n 



IUYl=v'2(rcf/p)i'''* 
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FIG. 2. Phase space: background counter-flow velocities U versus wave period T. The 
blocking curves U* are marked in thick black lines. The white-hole horizon (corresponding to 
the gravity-wave blocking curve Ug — thick red line — when 7 = 0) intersects the blue horizon 
or blocking curve for the blueshifted waves (the lowest thick black line, which asymptotes 
to Uy from below) and creates a cusp {Tc, Uc) below which both horizons disappear. Note 
that there also exists a blocking line for the negative-frequency waves, which intersects Ug 
at {Tb,Ub) and also asymptotes to Uj (from above). See also ^lOj. The four coloured stars 



correspond to the experimental results reported in Section IV 



to the appearance of a fictive symmetry in the representation space of the cubic-in-A; 
equation (pi). The resulting pitchfork bifurcation is then associated with the breaking of 
this symmetry, which is absent from the original system. In other words, below Tc, the 
blocking lines or horizons corresponding to both saddle-node bifurcations completely 
disappear. 

Note that the cusp is analogous to a critical point (second-order phase transition) 
in a thermodynamical phase diagram (Fig. |2]). The saddle-node lines (first-order phase 



9 

transitions) separate the analogues of thermodynamical phases pLQj. One can distin- 
guish seven regions of interest, marked by Roman numerals, which can be grouped into 
four phases: A = I + 11 + III, B = IV, C = V + VI and D = VII. 

The vertical line T = T^ roughly separates the capillary (T < T^) and the gravity 
(T > Tc) regimes. The A phase corresponds to a simple root of the cubic dispersion 
relation, where the incident wave can be of capillary (/) or gravity (// or ///) type. 
The saddle-node line ending on Uc (critical point) and Uj (tricritical point at infinity) 
corresponds to the threshold for the simultaneous appearance of blue-shifted waves and 
capillary waves propagating in the same direction as the incident ones {B). The C and 
D phases are characterized by the presence of negative energy waves. In the D phase, 
these negative energy waves are of the gravity and capillary type and propagate in the 
same direction as the counter-flow. The C phase is the forbidden region for gravity 
waves coming from A across B or D. Only gravity waves from the A phase (after mode 
conversion) or directly capillary waves are allowed to go into the C phase. 



III. EXPERIMENTAL SETUP 

We performed laboratory experiments to corroborate various aspects of the phase 
diagram in Fig. [2j These experiments were performed at ACRI, a private research 
company working on environmental fluid mechanics such as coastal engineering. The 
experiment features a wave-tank 30m long, ImSO large and ImSO deep, see Fig. [3] The 
piston-type wave-maker can generate waves with periods T =0.35-3s and typical wave 
heights of 0.5-30cm. A current can be created along or opposite to the direction of 
wave propagation with a maximum flow rate around 1.2m^/s. The waves themselves 
are recorded using several video cameras and the videos are digitalized and calibrated. 

In order to generate a gravity-wave horizon, a bump is immersed into the channel. 
The bump has a positive and a negative slope separated by a fiat section. We send 
a train of progressive water waves onto the bump, hindered by a reverse fluid flow 
produced by a pump. 

The background flow velocity depends on the water depth through flow rate con- 
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servation. The counter-current accelerates as the water height diminishes, reaches a 
maximum on the flat part of the bump before slowing down again as the bump height 
decreases. 




FIG. 3. Experimental setup. Characteristics of the bump (from left to right): linear slope 
with angle ai = 7.5° and length li = 8m, flat part I2 = 4.80m, linear slope 03 = —18.5° and 
I3 = 3.30m. Water depth (min-max): 50cm-160cm. 



IV. EXPERIMENTAL RESULTS 



We have performed detailed measurements at the {U, T) values marked by the four 
coloured stars in Fig. [2] Continuous low-amplitude wave-trains were used in order to 
minimize non-linear effects. The corresponding experimental space-time diagrams are 



11 



shown in Fig. |4| 

In the absence of a counter-current, the gravity and capillary terms in the dispersion 



relation ^ are equal for u = ^Jlgjlf. = 86rad.s~^ (/ = 13.7Hz), i.e. T = 0.073s. 
The range T > 0.1s corresponds (by convention) to a pure gravity regime, while a 
pure capillary regime exists for T < 0.04s. In the presence of a counter-current, one 
must look at the phase diagram in Fig. [2] to distinguish the gravity and capillary 
influence. The upper diagrams in Fig. |4] correspond to waves with a period of T = 
Is. These are therefore pure gravity waves, and we expect them to be blocked near 
\Ug\ = gT/SiT = 0.39m.s~^. The lower diagrams show the propagation of waves with a 
period T = 0.4s, approximately the lowest period allowed by the wave-maker. For weak 
counter-currents, these behave as pure gravity waves, but they should suffer a strong 
blueshifting towards the pure capillary regime as the counter-current increases, and 
penetrate through the gravity-wave white-hole horizon. The four diagrams in Fig. |4] 
thus correspond to the following cases. 

(a) First (red star), we recorded the normal propagation of a gravity wave of T = 
Is and amplitude A = 3cm against a moderate counter-current {\U\ increases from 
0.074m.s~^ to 0.087m.s~^ from left to right, well below \Ug\ = gT/S-rr = 0.39m.s~^), see 
Fig. g (top left). 

(b) Second (yellow star), still for gravity waves of T = Is and A = 3cm, we depict the 
range \U\ = 0.45-0. 55m. s~^ for which we recover the existence of a white hole marking 
a forbidden region into which gravity waves cannot enter (|6^g^''| ~ 0.53m.s~^), see 
Fig. Ill (top right), in agreement with the value measured in [221 12B]- We will come 
back to the apparent mismatch with the theoretical prediction |f/*^| = 0.39m.s~^ in 
the following section. Note that the white horizon or blocking line is actually blurred 
into a "blocking region". Also notice the clear blue-shifting due to the increasing 
counter-current: the slope of the world-lines increases from left to right, in full analogy 
with the behavior of light close to a gravitational fountain. Indeed, the slope is the 
inverse of the phase velocity and therefore proportional to the wavelength. The slope 
of incoming rays grows until the rays disappear at the horizon. We point out forcefully 
that the slope does not increase to infinity. This is due to the dispersive effect close to 
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the horizon, which leads to an Airy regularization mechanism, which we will discuss 
below. A "trans-Planckian" problem is thus avoided for the incident wave, since the 
blue-shifting does not become infinite, as in the purely relativistic case. However, as we 
will also discuss below (see Fig.lo]), there is still a problem for the mode-converted blue- 
shifted waves (as well as for the negative waves). In the pure gravity case, these would 
have fc — )■ cxD as f/ — )■ 0: if capillarity did not come into play, then the trans-Planckian 
problem would, in a sense, simply be displaced to flat space-time. 

(c) Third (blue star), we show the normal propagation of a capillary-gravity wave 
(T = 0.4s and A = 1 — 2cm) in the absence of a counter-current (f/ = 0), see Fig. |4] 
(bottom left). Note that there is an excellent agreement between theory and experiment 
about the wavelengths in this case: Aexp = 0.25m versus Ath = 0.249m in the whole 
range h = 0.50m — 1.60m. 

(d) Finally (green star), we show the propagation of a capillary-gravity wave (T = 
0.4s < Tc, A = 1cm) against a counter-current well above the blocking value: \U\ = 
0. 232-0. 275m.s~^ versus |f/g| = gT/8n = 0.156m.s~^, see Fig. (bottom right). From 
the diagram, we conclude that there is a complete penetration across this part of the 
white-hole region (forbidden for gravity waves), and barely any noticeable blue-shifting, 
contrarily to the case (b) of the blocked gravity waves. This last diagram therefore 
shows a double discrepancy with the theoretical expectation. 

First, the lack of any strong blueshifting with respect to the previous case means 
that the transition to the capillary regime has not been fully completed for these 
values of U. Indeed, for the counter-currents in the range of this camera position 
(|f/| = 0.233 - 0.275m.s-^), Ath = 8.95 - 6.13 x lO^^m whereas the measured Aexp is 
of the same order of magnitude as the incident wavelength. The capillary conversion 
seems to have taken place at a much higher value of the counter-current than expected. 
Indeed, further upstream (at higher values of the counter-current \U\), the appearance 
of capillary waves can be observed with the naked eye, see Fig. [5j We believe that the 
mismatch is due to the appearance of a transversal instability, see Fig. |6| which blurs 
the capillary conversion. We plan experiments with a narrower wave channel in the 
near future in order to reduce this transversal instability and study the conversion to 
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FIG. 4. Spatio-temporal diagrams for the four values of (U, T) marked with coloured stars in 
Fig. [2} The light and dark lines represent the world-lines of crests and troughs, respectively. 
The diagrams show: the normal propagation of a gravity wave against a moderate counter- 
current (top left); the blocking of a gravity wave at a blocking line or white-hole horizon due 
to a strong counter-current (top right); the normal propagation of a capillary-gravity wave 
in the absence of a counter-current (bottom left); and the propagation of a capillary-gravity 
wave across a region with a counter-flow velocity well above the gravity- wave blocking value 
(bottom right). The width of the images along the x-coordinate is 78.1cm (top left), 191.7cm 
(top right), 78.1cm (bottom left) and 113.92cm (bottom right). Other parameters: see main 
text. 
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the capillary regime in more detail. 

Second, in the absence of such a full transition to the capillary regime, the 
waves should have been blocked at (or near) the gravity-wave blocking velocity 
\Ug\ = qT/Stc = 0.156m.s~^. This blocking has not taken place either. This sec- 
ond discrepancy is in the line of the mismatch mentioned above in the pure gravity 
case, and we will come back to it in the next section. 




FIG. 5. Appearance of capillary waves for T = 0.4s in the presence of a strong counter-current 
U. 



V. AIRY INTERFERENCE AND GRAVITY- WAVE BLOCKING 



Airy interference provides the crucial mechanism through which a divergence of the 
amplitude is avoided [H] in the regime kh ^ 1, see Fig. [TJ An explicit expression for 
the stopping length Lg associated to the Airy interference in the case of pure gravity 
waves is 



Ls 



l9 



rp5/3 f[ 



(8) 



with x^^ the horizontal blocking position. A simple derivation of this expression can be 
found in the Appendix. Because of the geometry of our experiment, the background 
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FIG. 6. Development of transversal instability for T = 0.4s: near the wave-maker (at the 
furthest end of the channel), the wavefronts are nearly perfectly perpendicular to the channel's 
edges. Towards the bottom of the picture, the wavefronts start to deform under the effect 
of a transversal instability, and the free water surface acquires a 'fish-scale' pattern. The 
capillary waves in Fig. [5] appear on the fronts of these fish scales. 



surface velocity evolves linearly: |f/(a;)| = 0.51 — 0.05x, where x = corresponds to 
the kink in the mobile floor (note that the x-axis is oriented along the background 
flow, i.e. from right to left in Fig. Isl). From Fig. W\ the apparent mismatch between the 
well-known theoretical prediction Ug^ = — |^ for the counter- flow velocity at blocking 
{\Ug^\ = 0.39m.s"^ for T = Is), and the measured value (|f/g''''| = 0.53m.s~^) can also 
partially be understood. Ug^ is obtained in the ray approximation. The waves will 
actually be blocked a certain distance Ax^g^^^g further due to the Airy interference. 
The experimental blocking position xf^^ = —0.36m corresponds to the wave-blocking 
position. Taking the conservative assumption that a wave is considered "blocked" (i.e., 
it is no longer detected on camera) when its amplitude has decreased below 1% of its 
maximum, one obtains Ax^^^gg = 3.11Ls, i.e. (for T = Is) Ax^^^^gg = 0.61m. The 
ray-divergence position would then be at x^''^ + Ax^g^^^g = 0.25m, corresponding to 
a velocity f/^^P"'''^y = 0.50m.s~^. The theoretical blocking position corresponding to 
\Ug^\ = 0.39m.s~^ is x**^ = +2A2m. In other words. Airy interference explains roughly 
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20-25% of the difference between the theoretical prediction and the experimental mea- 
sure. 





rays 




imLs 



3.1138L, 



FIG. 7. Blocking of rays (in red) versus waves (in blue). In the ray approximation, the 
amplitude theoretically diverges towards the blocking point x* , beyond which it vanishes 
discontinuously. This divergence is regularized in the wave picture through Airy interference. 
The wave blocking point x^^ves li^s several Airy stopping lengths L^ further than x*ays. 



A second element which contributes to the difference between U^^ and U^^^ lies 
in the decrease with depth of the real velocity profile. The theoretical prediction for 
the value at the surface should, in a real experiment, be considered as an averaged 
(integrated) value over some depth, necessarily leading to a slightly higher value at the 
surface. The vertical velocity profile is approximately of the so-called plug type on the 
flat part of the bump and acquires a parabolic form after the flow has decelerated on 
the descending slope of the bump, see Fig. [8j Comparison of the surface values with the 
average, vertically integrated values gives differences of 5-20%. Interpolating between 
the cases represented in Fig. [SJ we obtain an estimated difference of ~ 10% between 
the theoretically predicted value and the value measured experimentally at the surface 
near the blocking point {xf^^ = —0.36m for T = Is). 

A third element which might be thought to be important is that the blocking velocity 
U^'^ = — 1^ is obtained in the pure gravity-wave limit, whereas even for T = Is, 
a small capillary influence persists and slightly increases the real blocking velocity. 
However, this difference is negligible, as can be seen from Fig. Ill where it corresponds 
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FIG. 8. Vertical background-flow velocity profiles: nearly plug-type profile on the bump at 
x=-2.43m (left) and parabolic profile towards the wave-maker at a;=4.20m (right), where x=0 
corresponds to the kink in the bump, see Fig. [3j 



to the departure between the red line Ug and the corresponding black line U* obtained 
numerically from the full Eq. ([s]). It can be estimated quantitatively as follows. In the 
regime klc 'ti 1, the dispersion relation ^ can be approximated by 



u 



Uk + 



gkh 



(kl. 



The condition ^ = for wave blocking becomes 



(9) 



U* 




|(i + i'?*^) 



(10) 



In the limit Ic = 0, the blocking wavenumber kg (corresponding to Ug = —g/Au) gives 



Kn 



Auj^ 



Inserting this value in the previous expression leads to 



U* 



U, + AU* = - 



3l_ 

Stt 



SOtt^ 



gT^ 



(11) 



For T = Is, At/* ~ -0.002m.s-^ < Ug. 

Non-linear effects could also play an important role, in spite of our attempts to limit 
them by working with small amplitudes. We limit ourselves to two comments. First, 
the finite wave amplitude A increases the blocking velocity. This well-known (but 



18 

poorly understood) phenomenon |22l El] can to a first approximation be modelled as 
an effective surface tension (see e.g. p9]): 

cu^ = gk{l + A^e) (12) 

(for pure gravity waves in the absence of a counter-current). Since A/lc ~ 10 in 
our experiments, it is clear that the finite amplitude has a much stronger infiuence 
on the blocking velocity than the intrinsic surface tension, and is perhaps the main 
contributor to the difference between U^^ and U^^^. This illustrates the importance of 
generating waves with a low factor Ak {Ak ~ 0.1 at T = Is in our current setup). A 
second non-linear effect which plays an important role in certain other wave-blocking 
experiments [231128] is the Benjamin-Feir instability, which leads to the appearance of 
so-called side-bands (excitations at frequencies slightly different from the fundamental 
one). However, we have verified conservation of period in our experiments, thereby 
nearly excluding this possibility. 

We believe the above elements to provide a reasonable explanation for the apparent 
mismatch between jf/g^'^l ~ 0.53m.s~^ and \Ug^\ = 0.39m.s~^. It should be noted that 
this mismatch is well known in the fiuid-mechanics community, but not well understood. 
For example, [22] and [28] also observe the blocking of T = Is waves at a countercurrent 
velocity \Ug^^\ ~ 0.53m.s~^, but do not attempt to interpret this discrepancy with the 
theoretical prediction. 

Also note that similar arguments would apply to the case of the gravity-capillary 
waves at T = 0.4s. There, however, the additional appearance of a transversal instabil- 
ity mentioned above further complicates matters and further experiments are required 
to clarify the situation. 

It is remarkable that the Airy interference hybridizes the character of the incoming 
wave. A hybrid is created between the original wave (through the period T, in the 
expression (pi) of the stopping length Ls) and the background fiow (through ^). In 
our experiment, we are sending continuous wave-trains. If one were to send wave- 
packets ("particles", i.e.: superpositions of waves), then these would be deformed into 
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superpositions of wave-flow hybrids, or "hybridons" . As a matter of fact, we can take 

1/3 

this observation further. The Airy stopping length obeys Ls oc X*^n , where A* is the 
wavelength at blocking, and ^„ the dimensionless "Unruh" number ^„ ~ ^ {^) 
obtained from the two characteristic "frequencies" at blocking: the wave frequency 
u and the flow gradient ^. This leads to the following interpretation. Dispersion 
has a double role in the near- horizon physics. It keeps the wavenumber finite (i.e., 
it solves the trans-Planckian problem — see next section), thereby avoiding the first 
relativistic ray-theory pathology A* — )■ 0. Second, it creates an interference mechanism 
which hybridizes this wavelength with the background flow by modulating it through 
^n- Dispersion thus replaces the wavelength by a characteristic interference length 
Lg. This mechanism of interferences allows to solve the second pathology associated 
to the ray theory: the infinite amplitude at the blocking point. Indeed, when ^n ^ 1, 
i.e. when the frequency of the wave is large compared to the spatial variation of the 
background flow velocity, then the WKB-approximation is valid. Note that ^„ ^ 1 
also leads to exp (cu/ (^) 1 — 1^1 and therefore to negligible Hawking radiation. 
Near the blocking point, though, one always has ^„ ~ 1. The WKB approximation 
then breaks down, and two resonance mechanisms come into play. The first one (Airy 
interference) is an adiabatic process and the second one (Hawking radiation) is a non- 
adiabatic process, see the Chapter on "The Basics of Water Waves Theory for Analogue 
Gravity" elsewhere in this Volume. 

As a final note, we should point out that we have neglected the presence of a zero 
mode in our considerations on the Airy mechanism. Such a zero mode (an u = 
solution to the dispersion relation, or superposition of various such solutions) would 
deform the free surface and complicate the interference pattern, see [HJ |30] and the 
discussion in the Chapter "The Cerenkov effect revisited: from swimming ducks to zero 
modes in gravitational analogs" elsewhere in this Volume. This omission is justified 
since in the regime kh ^ 1 one can minimize the amplitude of the zero mode by 
working at low velocities and limiting the slope of the bump. Note that 7 7^ implies 
the existence of a threshold \U\ > \U^\ for the appearance of a zero mode, contrarily 
to the pure gravity case. 



20 
VI. THE TRANS-PLANCKIAN PROBLEM 



Our experimental results to corroborate the theory developed in [TO] have been 
slightly marred by the appearance of a transversal instability, which we hope to rem- 
edy using a narrower wave-channel. Nevertheless, there is little doubt that capillary 
waves can penetrate through a gravity-wave blocking line. The full strength of this 
statement becomes clear in the context of the gravitational analogy. The surface ten- 
sion constitutes a high-fc dispersive correction to the low-A; gravity-wave theory. Such 
dispersive corrections can therefore completely alter the properties of a horizon: disper- 
sive horizons are no longer one-way membranes, and the infinite blueshifting associated 
with strictly relativistic horizons disappears. The example of surface waves shows that 
this statement can be true even if the dispersion is (initially) subluminal. 

The idea that dispersive corrections could solve the trans-Planckian problem of 
gravity has from the start been one of the cornerstones of the analogue gravity pro- 
gramme [6] . Most work has historically focused on the study of subluminal ( "normal" ) 
dispersion. This seems curious in the light of the following observations. In the pure 
gravity-wave (subluminal) regime, the trans-Planckian problem is indeed avoided for 
the incident wave: it is mode-converted into a blue-shifted wave and the wave bounces 
away from the horizon (in terms of the group velocity; the phase velocity is still directed 
towards the white hole). However, as this blue-shifted wave approaches "flat space- 
time" (i.e., as the counter-current velocity \U\ — )■ 0), a new, secondary trans-Planckian 
problem arises: the wavenumber of the blue-shifted wave ks ^ oo, see Fig. |9] This 
was observed earlier [31] with respect to Unruh's original subluminal model [32], and 
related problems with other subluminal models were also discussed in [33] and [M] . 
The same secondary trans-Planckian problem occurs for the negative-frequency waves 
associated to a Hawking-like process. Actually, even without invoking any horizon 
effects, a similar problem arises: Any counter-current, no matter how small, would 
allow for the existence of both blue-shifted and negative-energy waves with infinitely 
small wavelengths. Purely subluminal dispersion would then solve the primary trans- 
Planckian problem associated to the horizon, at the cost of creating a new one in flat 
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space-time. 




FIG. 9. Gravity waves and the trans-Planckian problem. For a given frequency to (dashed red 
Hne), a blue-shifted wave ks is created from the incident wave kj through mode conversion 
at the blocking line U = Ug. Since kB has Cg < 0, it moves away from the horizon towards 
lower \U\. As |C/| — )• 0, A;^ — )• oo, leading to a secondary trans-Planckian problem in flat 
space-time. 



More complicated dispersion relations, e.g. as in Helium-II (superfluid ^He), which 
has a "roton" minimum in the u'—vs—k diagram, could overcome this problem by 
letting the outgoing blue-shifted wave decay at the end of the quasiparticle spectrum 
into two rotons, with the same total energy and momentum (i.e., u — )■ uj/2 + u)/2; k — )■ 
k/2 + k/2) [31j. The secondary trans-Planckian problem is then avoided because the 
dispersion curve ends in such a two-roton decay channeljj 

Another obvious way of resolving the secondary trans-Planckian problem with sub- 

luminal dispersion is through dissipation, for example due to viscosity. Note that 

viscosity, apart from leading to dissipation, necessarily also introduces dispersion [37] . 

Dissipation in media is quite generic, often unavoidable (as in optical media, where dis- 

^ Actually, this is not the end of the story: these rotons are still subject to the background flow, and 
will therefore also start blucshifting, just like the original mode, and again split into two rotons each 
at the end of the quasiparticle dispersion curve, etcetera, leading to an apparently endless creation 
of rotons. This process is limited because the roton creation will deplete the superfluid component, 
and ultimately destabilizes the white-hole configuration. Also note that the rotons will relax after 
some time due to interaction with the environment and condense into a roton BEC |35j . If one 
takes such a ^He-like dispersion model seriously for true gravity, this could lead to the creation of 
a photon condensate near a white hole, or vice versa: outgoing particles from a black hole might 
originate from a condensate in curved spacetime. Although such a scenario might not be as crazy 
as it sounds [36] . we will not pursue this exotic line of thought further here, and stick to "simpler" 
solutions of the trans-Planckian problem. 



22 

sipation and dispersion are coupled through the Kramers-Kronig relations), and might 
be relevant for our "fundamental" space-time as well. Actually, from a theoretical QFT 
point of view, Lorentz symmetry violation is automatically accompanied by dissipation 
under quite general assumptions, and dissipative effects should therefore in principle be 
treated together with dispersive ones [M]- However, dissipation is formally much harder 
to treat than dispersion, and it has received little attention in the context of (analogue) 
gravity. Also, the relevance of dissipation obviously depends on its characteristic scale. 
Here, we mainly wish to stress that the secondary trans-Planckian problem in surface 
waves is solved through dispersion before dissipation becomes relevant. 

These observations suggest a more general interpretation for our results with water 
waves. The mesoscopic scale of surface tension "saves" the fluid continuum approxima- 
tion from breaking down in the presence of a counter-current: The capillary behaviour 
at high k is essential in order to avoid a trans-Planckian pathology. We therefore expect 
that any system displaying analogue gravity behaviour at low k through propagation 
on a moving background medium will necessarily have superluminal dispersion at suf- 
ficiently high fc's, unless dissipation kills the whole phenomenon before such scales are 
actually reached. The exotic case of Helium-II mentioned above might be the exception 
that confirms the rule, since even then, the dispersion relation is superluminal for a 
certain range of high wavenumbers beyond the roton minimum. 

Moreover, we can establish the following general rules. If the first corrections at 
intermediate k are subluminal, followed by superluminal corrections at high k, as in 
the case of deep-water waves, then there will always be a white horizon and a blue 
horizon, leading to the possibility of a Badulin-type double-bouncing scenario. There 
will then also always exist some counter-current velocity Uc for which these white and 
blue saddle-node points merge into a pitchfork bifurcation and both horizons disap- 
pear, allowing for direct dispersive penetration. Such direct dispersive penetration is 
actually even more universal: it suffices to have a superluminal correction at high fc to a 
relativistic low-fc behaviour, irrespective of the intermediate regime. This is true, e.g., 
for phonons in BECs [39], just like for capillary surface waves: sufficiently blueshifted 
( "superblueshifted" ) modes will always be able to penetrate through any counter- flow 
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barrier, unless dissipation prevents such superblueshifting. The case of surface waves 
is in a sense richer than that of BEC-phonons, in that there is a true blocking line for 
low-frequency gravity waves, which cannot directly penetrate the horizon. For BECs, 
the absence of an intermediate subluminal correction implies that even low-frequency 
phonons can in principle superblueshift and cross the horizon directly. Finally, in spite 
of our several comments regarding dissipation, it seems that both horizon-crossing 
scenarios can indeed be fully realized for surface waves before being dissipated. 

To sum up, to enter a white hole — or, by time-inversion: to escape a black hole — one 
has to either tune the period to be subcritical, or bounce on two horizons. 

The bottom line is of course what this teaches us for real gravity. Here the issue 
is more complicated, because extrapolation from the current state of observations has 
so far not given any evidence for dispersion (or dissipation) even at the Planck scale. 
Other complications might also arise which are peculiar to real gravity. For example, 
in |l0] it was shown using a simple toy model that superluminal dispersion would render 
gravitational black holes strongly unstable, due to the leaking of resonant modes. In 
any case, we may conclude that, if dispersion is indeed relevant in gravity (possibly even 
far beyond the Planck scale), then subluminal dispersion alone would most certainly 
not be sufficient to solve the trans-Planckian problem. 



APPENDIX: AIRY STOPPING LENGTH 

Smith was the first to derive the Airy equation in 1975 by performing an asymp- 
totic expansion of the water-wave equations (Euler equations + continuity equation 
+ boundary conditions) close to the caustic [19]. He inferred the so-called amplitude 
equation which is a nonlinear Schrodinger equation with a term proportional to the 
distance to the caustic. When the cubic term is negligible, the amplitude equation 
reduces to the Airy equation. 

Following Smith and after tedious algebra (matched asymptotics and WKB solu- 
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tions), Trulsen and Mei computed the following stopping length (for 7 = 0) in 1993 [41]: 

In 1977, Basovich & Talanov |21] provided another derivation of the Airy equation 
by noticing that ^ = at the blocking point. Taylor-expanding the function U{k) 
close to its parabolic minimum and the function U{x) close to the stopping point x=k 
and taking the inverse Fourier transformation, they deduced the Airy function and the 

associated stopping length: 

1/3 

" * (14) 



In 1979, Peregrine & Smith [20j used an operator expansion method: the idea is to 
inverse Fourier-transform a truncated series expansion of the dispersion relation written 
in the form G{uj^ k,x) = 0. Again, the cubic Schrodinger equation with a spatial term 
was derived with another expression for the stopping length: 



2G, 

where the subscripts mean partial derivative and the derivatives are taken at the block- 
ing line. The method was generalized in 2004 by Suastika PH [25] to include viscous 
dissipation and wave breaking. In 2003, Lavrenov |17j applied a saddle-point method to 
the Maslov integral representation of the uniform wave field asymptotics in the vicinity 
of the blocking line. He found: 

"fcfe 



^x 



^-lifl (^«) 



where u = Q{k,x) is the dispersion relation function and the derivatives are taken at 
the caustic. 

We will show that it is possible to derive the stopping length in a very simple fashion, 
inspired by the method of Basovich & Talanov, and derive a previously unnoticed 
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scaling law for Lg. 

We write the background flow velocity near the critical value ?7=k = — |^ as a function 
of X and /c, and develop both to lowest non-zero order around the stopping length: 

U{k) ^ U{K) + ^^^(A; - Kf = U,- ^{k - k,f (17) 



U{x) ~ t/, + ( ^ ) (x - X,) (1^ 



2 ' ^ 4^2 

Equating both into 



and making the substitution H{x) ~ Qiik-k*)x jniniediately leads to the Airy differential 
equation 

g-A-i/ = (20) 

where X = ^r^, with L, = '— — t-r- Thus, H(x) is an Airy function 

i7(^) ^ A, (^^) = - /" cos T-t^ + ^^t ) rft (21) 

and 

-1/3 



^--(f)::. 



' 16(27r5)i/3- 

is the Airy stopping length, which depends both on the incident wave and the back- 
ground flow: it scales with the period T of the incident wave as Lg oc T^/^ and with 

-I /O 

the background flow acceleration as L^ oc (^) _ . A straightforward dimensional 
analysis {Lg ~ g"T^ (^) _ ) would only lead to a = 1 and /3 — 7 = 2. 

Note that Airy interference requires the flow gradient to remain approximately con- 
stant over the characteristic length of the interference process. This can easily be seen 



in our derivation: Eq. (18) is only a good approximation if ^-^ f« for x — x^: = 0{Ls 



In our experiments, ^ is a constant on the whole linear slope of the bump where the 
horizon x* is located, so we do not need to worry about this issue. 
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